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SUMMARY 


An analytical design procedure for Leading-Edge Extensions (LEE) 
has been developed for thick delta wings. This LEE device is designed 
to be mounted to a wing along the pseudo-stagnation stream surface 
associated with the attached flow design lift coefficient of greater 
than zero. The intended purpose of the device is to improve the 
aerodynamic performance of high subsonic and low supersonic aircraft at 
incidences above that of attached flow design lift coefficient, by using 
a vortex system emanating along the leading edges of the device. The 
low pressure associated with these vortices would act on the LEE upper 
surface and the forward facing area of the wing leading edges, providing 
an additional lift and effective leading edge thrust recovery. 

The first application of this technique was to a thick, round- 
edged, twisted and cambered wing of approximately triangular planform 
having a sweep of 58° and aspect ratio of 2.3. The panel aerodynamics 
and vortex lattice method with suction analogy computer codes were 
employed to determine the pseudo-stagnation stream surface and an 





optimized LEE plan-form shape, respectively. 

The aerodynamic effectiveness of thirty six different LEE planform 
shapes were examined for the given wing by considering the influence of 
geometrical parameters such as chord, sweep angle and span extent. This 
investigation showed that the outboard reduction of the LEE span-extent 
minimizes the lift-to-drag ratio, regardless of the LEEs' planform-shape 
and area. Also, with the same planform area, it was found that constant 
chord is relatively more effective than LEEs having sweep angles less 
than that of the wing. Further, relative to the wing root chord and 
span, a 3.2% constant chord LEE with 89% span extent was selected as 
being the best candidate for the final design planform. 
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Chapter 1 


INTRODUCTION 

Future high-subsonic or supersonic cruise swept wing aircraft are 

likely to be required to operate efficiently over an extended portion of 

their flight envelope. There are two basic approaches for designing 

such aircraft. The first, is a conventional approach, and seeks to 

maintain fully attached flow at each point of the envelope, whereas, the 

second approach attempts to use the organized separated flow at off- 

design and attached flow at design conditions. The design criterion of 

the conventional approach is more desirable, because an aerodynamical ly 

efficient aircraft always achieves its best performance with attached 

flow unless the wing is extremely slender. The primary cause of this 

high efficiency is the production of aerodynamic thrust associated with 

attached flow at the leading edge, as well as in the warped camber. The 

schematic flow representation of the techniques used to retain attached 

flow conditions at the leading edges of such swept wing aircraft are 

shown in Fig. 1.1. Variable camber at the leading edge, leading-edge 

flap, and large leading-edge radii are known for their potential to 

* 

delay the onset of the leading-edge flow separation [1-3]. However, the 
♦The numbers in brackets indicate references. 
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natural tendency of flow towards separation for these wings, especially 
at off-design conditions such as for take-off, landing, and maneuvering, 
appears inevitable (see dash lines in Fig. 1.1). At off-design perform- 
ance the flow characteristics of such aircraft are changed dramatically 
by the formation of a generally stable and coherent leading-edge vortex 
system. These vortices, which result from the leading-edge flow sepa- 
ration and subsequent flow reattachment on the remainder of the wing, 
are responsible for the changes which occur in the aircraft aerodynamic 
characteristics. For the purpose of illustration. Fig. 1.2 shows the 
typical flow types occurring over such aircraft at design and off-design 
conditions. 

The resultant vortex system generates additional lift, caused by 
the low pressure regions under the stable vortex system, and produces 
the well known nonlinear aerodynamic behavior called "vortex lift." A 
typical comparison of the vortex lift, to that of attached flow lift is 
shown in Fig. 1.3. Accompanying the additional lift is the increased 
drag which results from the loss of the leading-edge suction associated 
with attached flow around the leading edge [4]. In addition to the .drag 
penalty, the inboard movement of the center of the vortex coupled with 
the flow failing to reattach with increasing angle of attack, and vortex 
breakdown, results in a pitch-up condition [5]. These characteristics 
restrict high-g subsonic and transonic sustained maneuver, because of 
the excess engine thrust needed to overcome the drag-due-to-lift. As a 
result, the aerodynamic characteristics of the naturally occurring 
leading-edge flow separation over swept wing aircraft operating at off- 



a. variable camber b. leading edge flap c. blunt leading edge 


Fig. 1.1 Some techniques used at the leading edges of an aircraft to attain 
attached flow. (Streamwise cut) 



Attached flow Vortex flow 

(Design condition) (Off-design condition) 


Fig. 1.2 Typical flow types occurring at the wing leading edges of 
an aircraft. 



Fig. 1.3 Lift characteristics associated with attached and separated 
flow. 
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design conditions must be considered early in the aircraft design 
cycle. 

As technology in aircraft design has developed, methods for improv- 
ing multimission capability have been explored. One such method, the 
subject of this study, is to design the wing to achieve fully attached 
flow at the cruise design condition, and controlled leading-edge sepa- 
ration at take-off, landing, and maneuvering [3]. This method is an 
alternative approach to the conventional attached flow, for designing a 
high-subsonic or supersonic cruise swept wing aircraft. The basic con- 
cept of this approach is to let the flow separate and roll up into an 
organized leading-edge vortex system which is located appropriately. 

For this purpose, a family of vortex control devices such as fixed and 
movable leading-edge extensions has been developed. For example, the 
Sharp Leading-Edge Extension (SLEE) is a fixed leading edge extension 
composed of a flat or bent plate attached to the wing lower surface and 
been employed on swept wing models with round leading edges [6, 7]. An 
example of a movable leading-edge extension is the Leading-Edge Vortex 
Flap (LEVF). This leading-edge device can be rotated about its hinge- 
line and set at scheduled deflection angles which vary with angle of, 
attack and Mach number [8, 9]. The flow mechanism associated with the 
leading edge for swept wing aircraft with SLEE and LEVF is shown sche- 
matically in Fig. 1.4. Unlike the conventional attached flow approach, 
SLEE and LEVF benefit from the natural tendency of flow separation at 
the leading edges. Such devices, when properly designed and positioned, 
can confine the entire leading-edge vortices to their upper surface and 



provide flow reattachment on the wing along the device knee or hinge 
line. As a result, the aircraft not only produces additional lift, but 
it also generates a thrust force component, as the low pressure associ- 
ated with the confined vortices acts on the neighboring surfaces. These 
leading-edge vortex control devices have been validated experimentally 
through extensive parametric studies on different wing models [5-10]. 

Following the latter approach (i.e., letting the flow separation 
occur), the present study attempts to develop a different expression of 
the fixed leading-edge extension device concept. The intended purpose 
of the device, designated as Leading-Edge Extension (LEE), is to improve 
the aerodynamic performance of high-subsonic and low supersonic aircraft 
away from the conditions for cruise in general called off-design. The 
scope of the present study will be discussed subsequently in this 
chapter. 

Pursuing the concept of leading-edge vortex control, a literature 
survey was conducted for devices that had potential for controlling the 
leading-edge flow separation of wings with moderate-to-highly sweptback 
leading edges. The next section presents an overview discussion on the 
aerodynamic effectiveness of devices such as SLEE and LEVF. 

1.1 Literature Survey 

In recent years, leading-edge extension devices have been the 
subject of extensive studies for improving the aerodynamic and pitching 
moment characteristics of high-subsonic and transonic aircraft, capable 
of maneuvering efficiently at high lift. One of the earliest efforts 
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made in using a leading-edge device was by Wilson and Lovell [10]. The 

objective of their study was to increase the C. of a 15% thick, 

L,max ’ 

blunt leading edge of an approximately 60° delta wing, designated as DM- 

1, by producing vortex flow over its upper surface. They determined 

experimentally that the flow separation required could be best achieved 

by employing a fixed part-span leading-edge extension on the vehicle. 

(The vehicle had a symmetrical airfoil and no twist so the design lift 

coefficient (C L d ) was zero). As a result, their wind tunnel studies 

indicated that, attachment of this device increased the oriqinal C, 

L,max 

of the DM-1 from 0.6 to 1.01 with essentially no drag penalty at low 
lift coefficients, only a slight increase at moderate values, and a 
significant drag reduction at lift coefficient beyond 0.75. This 
initial study enlightened the importance of a leading-edge flow control 
device in the aerodynamic performance of an aircraft. 

A series of different SLEE and LEVF devices have been investigated 
by Johnson and Rao [6] and Tingas and Rao [7] on a 60° swept, cropped- 
delta wing with round leading-edges. These experimental studies have 
shown a substantial improvement in drag reduction potential of both , 
devices at moderate-to-high angles of attack. These authors have con- 
cluded that there is a need for a concerted effort to optimize the 
effectiveness of these devices, which function by maintaining the vortex 
on the upper surface of device, flow reattachment where the device joins 
the wing and attached flow on the rest of the wing. These studies 
further identify the following methods as having significant potential 
to increase the aerodynamic efficiency of such devices. They are: 



tapering, twisting, segmenting, and proper mounting position on the 
wi ng . 

Additional details concerning leading-edge flow control devices can 
be obtained from recent publications by Lamar and Campbell [11], and Rao 
[ 12 ]. 


1.2 Present Study 

The objective of this study is to develop an extension to the 
device used by Wilson and Lovell, which would improve the aerodynamic 
performance and pitching moment characteristics of cambered and twisted 
high-subsonic and low-supersonic aircraft at off-design conditions. 

This leading-edge device, designated as a Leading-Edge Extension (LEE), 
is to be mounted to a wing along the Pseudo-Stagnation Stream Surface 
(PSSS) associated with the attached flow design lift coefficient (C^ j) 
of greater than zero [11]. The PSSS is a dividing stream surface which 
separates the incoming flow into two regimes, in general over the upper 
and under the lower wing surface. Two streamwise cuts through the PSSS 
are shown schematically in Fig. 1.5 to illustrate the surface curvature. 
The present study seeks to determine a representation of the PSSS based 
on the following four assumptions. 1) There exists a PSSS associated 
with a swept wing aircraft at attached flow design condition. 2) The 
intersection of the PSSS with a number of parallel xz planes spanning 
the wing produces curves. 3) These curves can be represented by the 
pseudo-stagnation streamline leading to the pseudo-stagnation points and 
derived from the local slopes of the resultant velocities (Vv| + »|) 
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Fig. 1.4 Schematic mechanism of leading-edge flows over a wing with 
SLEE and LEVF. (Streamwise cut) 


Pseudo-stagnation 
Streamlines ' 



I v xl H v zl 
|Vy| t 0 


Fig. 1.5 Schematic representation of the PSSS corresponding to two 
airfoil sections of a wing. 


at appropriate points in the xz plane. 4) A spanwise surface fitted 
linearly through the resulting intersections is an approximation of the 
PSSS described in assunption (1). 

In order to accomplish the task of designing an aerodynamical ly 
efficient LEE planform shape, an analytical procedure had to be develop- 
ed. This procedure which forms the basis of the present study can be 
outlined into two major steps: 

a) Analytical determination of the PSSS at attached flow design 
condition for the wing. 

b) Analytical optimization of the chordwise extent and the plan- 
form shape of the PSSS at separated flow conditions. This step 
would in fact determine the optimum LEE size for the given 
wing. 

To demonstrate the procedure outlined above, a candidate wing and 
computer codes (i.e., analytical tools) had to be selected. As a 
result, a thick, round-edged, twisted and cambered wing of approximate- 
ly triangular planform having a sweep of 58° and aspect ratio of 2.3, 
was chosen to provide the first application of this technique. The 
planform view of this wing model is shown in Fig. 1.6. At the outset, 
four computer codes were considered to be the options for analytical 
execution of the present study at a high subsonic Mach nunber. These 
codes were Free Vortex Sheet (FVS), Panel Aerodynamics (PAN AIR), Vortex 
Lattice Method with Suction Analogy (VLM-SA) [13-15], and a transonic 
computer code. Although, attempts were made to obtain and employ a 
transonic computer code in this study, due to the high subsonic Mach 
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numbers of interest, none was available when this study began which 
could reliably estimate the pressures on thick-delta wings. Following 
elimination of this code, the FVS was also excluded from the code 
options because of the authors' unsuccessful past experience which 
included efforts to obtain a converged solution for the DM-1 + LEE com- 
bination of Reference 10. Hence, the available computer code options 
were reduced to two, namely PAN AIR and VLM-SA codes. Further, after 
unsuccessful attempts for determining the velocity field solutions by 
VLM-SA because of the thickness omission by the code, the PAN AIR code 
was assigned to perform the task. The task was successfully accomplish- 
ed and the velocity field solutions for different wing sections were 
analytically determined at the attached flow design lift coefficient 
(C[_ d ) of 0.25 and Mach number of 0.8. (Note that this value of C|_ ^ 
was used to simulate the design angle of attack (a d ) of 6.0°). Neglect- 
ing the sidewash (V y ) effect, the resultant velocity vectors obtained 
from vectorial addition of the axial (V ) and the upwash (V ) velocity 
components associated with each wing section were plotted, and the 
corresponding pseudo-stagnation streamlines were graphically determined. 
The resulting stagnation streamlines were designated as being "pseudo" 
because they did not correspond to the actual stagnation point where the 
magnitude of the three velocity components are all zero. In fact, 
except at the center line of a three-dimensional swept wing there exists 
no other point on the wing surface, from a potential flow viewpoint, 
where zero sidewash velocity will occur. As a result, the pseudo- 
stagnation streamline solutions correspond to certain points on the wing 


surface where only the magnitude of V x and V z components of the total 
velocity are zero. Consequently, the determined pseudo-stagnation 
streamlines, are the planar cuts through a surface which represent the 
"pseudo" stagnation stream surface. Lastly, a portion of the determined 
PSSS solution is to be designated as the shape of the LEE device. 

The LEE acts as a dividing stream surface. In general, depending 
upon the accuracy with which the LEE (or the PSSS) is determined, its 
presence, other than a small skin friction drag, should not affect the 
main wing aerodynamic performance at the designed angle of attack (a^) 
of 6.0° (see Fig. 1.7a). However, at higher angles of incidence, vor- 
tices would be generated as a result of forced flow separation by the 
sharp leading edges of the LEE device. These vortices can be controlled 
through LEE planform shape optimization by varying parameters such as 
the chordwise extension, spanwise extension, and leading-edge sweep 
angle. As shown in Fig. 1.7b, a properly designed LEE planform, can 
capture the entire leading-edge vortices on its upper surface and 
provide flow reattachment at, or near the wing upper surface leading 
edge. The confined leading-edge vortex system induces suction pressure 
which acts on the LEE surface and the forward-facing area of the wing 
leading edges, providing an additional lift and an effective leading- 
edge thrust recovery. As a result, the aerodynamic thrust force gener- 
ated in the flight direction yields a reduction in drag, relative to a 
planar configuration, and the added lift permits the aircraft to operate 
at lower angles of attack which may delay the pitch-up moment problem. 
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Fig. 1.6 Planform of the wing model. 




Fig. 1.7 Leading edge flow mechanism of a wing-LEE combination at 
a) design angle of attack, b) off-design angle of attack. 
(Streamwise cut) 


Chapter 2 


DETERMINATION OF PSEUDO-STAGNATION STREAM SURFACE (PSSS) 

This chapter discusses in detail the application of the PAN AIR 
code and the procedure employed for determining the PSSS of the wing 
model at the attached flow design condition. Further, Appendix A 
provides a brief discussion on the PAN AIR code and a test which was 
conducted for validating the resulting PSSS solution. 

2.1 Method Employed 

PAN AIR code is a system of computer programs for the detailed 
analysis and the non-iterative design of arbitrary aircraft configura- 
tions in three-dimensional, steady, inviscid, irrotational, subsonic and 
supersonic flows. The configuration surface is partitioned into several 
"networks," each approximated by a set of panels on which unknown source 
and doublet singul arity distributions are assigned. By imposing bound- 
ary conditions at a discrete set of points, the integral equation 
solution to the partial differential equation is reduced to a system of 
linear algebraic equations relating the unknown singularity strengths 
which in turn determine the properties of the flow field. Additional 
discussion on some of the capabilities and limitations of the PAN AIR 
code are explored in Appendix A.l. 

The PAN AIR computer code was employed in the present study to 
determine the velocity field solution of the wing model at the attached 
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flow design angle of attack (a d ) of 6.0° and Mach nunber of 0.8. For 
this purpose, the survey network coupled with the Influence Coefficient 
(IC) update capability (see Appendix A.l) of the code was first exercis- 
ed. Several attempts were made, but due to the problems introduced by 
the IC-update package, efforts to use this economically efficient 
approach had to be terminated. As a result, only the survey network 
capability of the code was employed in the present study to determine 
the velocity field solution of the wing model at the design condition. 
This task was successfully accomplished using the PAN AIR code and 
resulted in graphical determination of the PSSS for the wing model. 
Sections 2.2 and 2.3 discuss the procedure involved in determining the 
PSSS. 


2.2 Survey Networks 

The survey networks adopted in the present study were vertical xz 
planes located at sixteen different stations along the semi-span of the 
wing model. These survey networks were generated such that, each would 
enclose the nose portion of its corresponding station and stand-off from 
the section a distance of approximately equal to .08% of the wing c r . 

The networks began at the upper surface just behind the leading edge and 
extended around the nose to the lower surface mid-chord. Due to the 
similarity of the survey network geometries and the involved process of 
their generation, only a typical survey network (located at the fourth 
station) shown in Fig. 2.1, will be discussed. This figure also shows 
the planform distribution of the other survey network locations over the 



semi-span of the wing model. Further, the enlarged cross sectional view 
of the survey network and the nose portion of its corresponding airfoil 
section at the fourth station is shown in Fig. 2.2a. Since the PAN AIR 
Code velocity field solutions were assigned to be calculated at the 
center point of each panel in a particular survey network, it was 
essential to provide the survey networks with enough panels, so that, 
once the resultant velocity vectors, associated with V x and V z , were 
plotted, the pseudo-stagnation streamlines could be depicted graphically 
for each wing section. For this purpose, a geometrical computer code, 
called GEOMABS [16], was employed to intensify the paneling on the 

survey networks. Figure 2.2b, shows the repaneled survey network. It 

can be seen from the figure that the panel density is concentrated 

primarily around the portion of the survey network which faces the nose 

of the associated wing section. This would provide more velocity vector 
solutions needed to determine graphically the accurate location of the 
resulting pseudo-stagnation streamlines, as they meet their 
corresponding wing section. Similar survey networks were generated for 
all sixteen semi-span stations of the wing model. Each individual 
survey network was positioned on the wing model and a separate PAN AIR 
code execution was performed. 

2.3 Flow Field and the PSSS Solution 
The PAN AIR code executions yielded the axial (V x ), the sidewash 
(V ), and the upwash (V ) velocity components at the panels center 

y 

points of each survey network. The resultant velocity vectors, obtained 
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Fig. 2.1 Semi-span planform view of the wing model with a survey 
network located at fourth station. 


mm 


Fig. 2.2 Enlarged cross-sectional view at fourth stati 

a) original panels. 

b) densed-up panels. 








from the vectorial addition of the axial and the upwash velocity compo- 
nents were plotted at the center point of each panel for a given survey 
network. An example of these plots is shown in Fig. 2.3a. For a given 
survey network, the streamline associated with minimum velocity magni- 
tude (i.e., |V |»|V |»0, pseudo-stagnation point) was drawn tangent to 
the plotted velocity vectors. As mentioned earlier, these streamlines 
were designated as pseudo-stagnation streamlines. Fig. 2.3b shows the 
nose portion of an airfoil section with its corresponding velocity field 
and the graphical pseudo-stagnation streamline solution. These graphi- 
cal streamline solutions yielded their coordinate points relative to the 
corresponding wing section. A cubic spline curve was fitted through the 
graphically determined coordinate points, associated with each pseudo- 
stagnation streamline, to ensure the smoothness of the resulting stream- 
line solutions. Each of these solutions were equally extended out a 
distance of 4.8" (i.e., 19% of the wing c f ) ahead of the wing leading 
edge. This distance was thought to be sufficient for the present analy- 
sis of the LEE device. 

The unrealistic velocity field solutions obtained at the tip region 
prevented the graphical generation of the pseudo-stagnation streamlines 
for the last two wing sections because of the manner in which the tip 
thickness was modeled.^ As a result, this unrealistic solution at the 
tip section affected the flow field of the neighboring station as well 

^It was subsequently learned that unrealistic flow field solutions could 
be expected to occur in the tip region, if the wing thickness there was 
zero (private communication with Larry L. Erickson of NASA Ames). It is 
anticipated that this problem will be reexamined in the near future, so 
that the remaining pseudo-stagnation streamline solutions, at the tip 
region, would be determined. 




tip section affected the flow field of the neighboring station as well 
(i .e., 15th station) . 

A warped surface was linearly fitted spanwise through the available 
pseudo-stagnation streamline solutions and was designated as the PSSS. 
The three views of the determined PSSS solution are shown in Fig. 2.4. 
Further, five sectional cuts through the wing-PSSS combination and the 
enlarged cross sectional view of the same cuts are shown in Fig. 2.5. 
Approximately the resulting PSSS has a semi-span of 14.33" (i.e., 89% of 
the wing semi -span) and a constant chord of 4.8". 

It was essential to examine the degree of accuracy of the deter- 
mined PSSS solution. For this purpose, the PAN AIR code was employed 
once again to model the wing-PSSS combination at the design condition by 
specifying the PSSS as being a lifting surface. This investigation is 
discussed in detail in Appendix A. 2. From this pressure distribution 
study, it was concluded that the resulting PSSS solution is a good 
approximation of the actual dividing stream surface associated with the 
wing model at the attached flow design condition. 
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Fig. 2.5 Sectional cuts through the wing-PSSS combination and enlarged 
cross sectional views. 
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Chapter 3 


RESULTING LEADING-EOGE EXTENSION (LEE) EFFECTS 

This chapter discusses in detail the application of the VLM-SA code 
in the present study for determining an effective LEE planform shape for 
the wing model. A brief discussion on some of the code capabilities and 
limitations is presented in Sec. 3.1, and these have been further 
expanded in Appendix B. The effect of various LEE constant chord, span, 
and sweep angle are discussed in Secs. 3. 2-3. 4, respectively. 

3.1 Method Employed 

Vortex Lattice Method coupled with Suction Analogy (VLM-SA) 
developed at NASA Langley Research Center estimates overall forces and 
moments of complex planforms at subsonic speeds. The code is based on 
steady, inviscid, irrotational , incompressible flows, and uses the 
Prandtl-Gl auert rule to account for compressibility. It approximates 
the continuous distribution of bound vortices over a lifting surface by 
a finite number of elemental panels which are then replaced by horseshoe 
vortices. The resultant force contribution of an individual panel is 
determined by imposing the no-flow penetration boundary condition to 
each of the elemental panels. These forces are then summed appropriate- 
ly to obtain lift, drag, thrust and pitching moment. Additional discus- 
sion on VLM-SA code is presented in Appendix B.l. 

The VLM-SA code is employed in the present study to investigate the 



effects of the presence of the LEE device, as well as its geometrical 
parameter variation, on the aerodynamic characteristics of the thick, 
twisted and cambered, basic wing model. Although the twist and camber 
of the wing model is represented by its mean camber surface, the thick- 
ness effect is ignored by the VLM-SA code (see Appendix B.l). A comput- 
er program was developed to generate the required slopes at the control 
point (also called the local angle of attack) of each elemental panel 
located along the mean camber surface of the wing model. This program 
was further modified and used to find the local angles of attack for the 
warped surface of the LEE device. These two programs are listed in Ap- 
pendix B.2. In addition, an effort has been made to evaluate the cap- 
ability of the VLM-SA code in predicting the total aerodynamic vortex- 
induced forces for a wing-LEE configuration, similar to the one employed 
in the present study. This discussion is presented in Appendix B.3. 

The analytical solutions for the basic wing model of the present 
study were first intended to provide a base line for comparative assess- 
ments of the LEE device. However, as discussed in Appendix B.3, this 
approach appears unjustified, because it is suspected that the VLM-SA 
would estimate the drag pessimistically in the case of wing-LEE analy- 
sis of the present study. As a result, throughout this study the aero- 
dynamic effectiveness of different wing-LEE combinations will be 
emphasized relative to one another rather than to the basic wing model. 

3.2 Constant Chord 

Analytical estimates of the total vortex-induced forces and moments 
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were generated on the wing model as well as with selected constant chord 
LEE configurations by employing the VLM-SA code. There were a total of 
six constant chord LEEs examined in this study. The selected chord di- 
mensions included 4.8", 3.6", 2.4", 1.6", 1.2", and 0.8". All these 
examined LEEs had a semi-span of 14.33". The planform view of the wing 
model with these LEEs are schematically shown in Fig. 3.1. The code 
estimates were generated using seven chordwise and 25 spanwise horseshoe 
vortices on the half-span of the wing model, and another 7-by-21 array 
of horseshoe vortices were used on the LEE planforms. 

The VLM-SA code estimates of drag and lift coefficients for the 
basic wing as well as for the wing-LEE configurations are plotted in a 
drag polar form in Fig. 3.2. Although, the solutions include the aero- 
dynamic performance of wing-LEE combinations below d (i.e., 0.25), it 
is practical only to examine their aerodynamic effectiveness at higher 
lift coefficients. As shown in the figure, at the design lift coeffi- 
cient, the attachment of the LEE device produces only a slight addition- 
al drag as compared to the basic wing, regardless of the LEE planform 
area. As was verified by PAN AIR code (Appendix A. 2), these results 
also indicate that the LEE (PSSS) devices are good representations of 
the dividing stream surface associated with the wing model at the 
attached flow C L d . These results are remarkable considering the dif- 
ferences between the theoretical methods employed by each code. 

The drag polar comparison reveals that 0.8" and 1.2" constant chord 
LEEs produce nearly the same aerodynamic effectiveness through lift co- 
efficient range of about 0.25 to 1.20. However, it appears that consid- 
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Fig. 3.2 Effect of constant chord length for LEE on drag polar 
niEE= 89%. 





erable improvement can be achieved in the lift and drag characteristics 
of the wing-LEE combination by employing a longer LEE chord extension. 
For example, as compared to 1.2" and 0.8", 2.4" constant chord LEE pro- 
duces less drag in the lift coefficient range of 0.50 to 1.20. Also, 
this figure shows for the same lift coefficient range, 4.8" constant 
chord LEE has a remarkable drag reduction capability. In fact, it pro- 
duces minimum drag beyond the lift coefficient of 1.0. Further, the 
same results are shown in Fig. 3.3, where lift-to-drag ratio is plotted 
against lift coefficient. This figure also indicates that, a 4.8" 
constant chord LEE achieves the best aerodynamic performance throughout 
the lift coefficients beyond ^ However, from the practical point of 
view it should be mentioned that the final LEE planform design should 
have a chord dimension which is relatively shorter than the wing local 
chord, especially in the tip region. A smaller chord LEE not only 
benefits from the reduced structural weight, but it also minimizes the 
effect of bending moment about the wing-LEE junction. This bending 
moment occurs at off-design performances as the low pressure associated 
with the leading-edge vortices act on the upper surface of the LEE 
device. As a result, the present study seeks to design a LEE which 
employs both minimum area and chord. These two design criterion are 
referred to as the design requirements of the present study. 

3.3 Sweep Angle 

The aerodynamic effects of six different LEE sweep angles (A^) 
were investigated in the present study by employing the VLM-SA code. 
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Fig. 3.3 Effect of constant chord length for LEE on lift-to-drag 


A schematic planform view of the wing model and the LEE sweep angles are 
shown in Fig. 3.4. The selected angles include 53°, 54°, 55®, 56®, 57®, 
and 58®. These angles are measured from a horizontal line which passes 
through a point located at a distance of 0.8" ahead of wing leading edge 
along the pseudo-stagnation streamline associated with the third semi- 
span station. All the selected LEEs had a semispan of 14.33" which is 
equal to 89% that of the wing model. Figure 3.4 also shows the LEE's 
tip chord dimensions. The VLM-SA estimates were obtained by using the 
same nimber of horseshoe vortices on the wing and the LEE planforms as 
were used in the previous section. 

The resulting VLM-SA solutions of lift and drag coefficients are 
shown in Fig. 3.5. This figure shows that, except in the lift coef- 
ficient range of about 0.6 to 1.0, variation has only a slight 

effect on the drag and lift characteristics of the wing-LEE configura- 
tion. This is an interesting result. As shown in Fig. 3.4, the LEE 
area and the corresponding tip chord decreases as A^ increases. In 
fact, 53® sweep angle LEE has twice the area as 57® sweep angle, 
however, they both produce almost the same drag characteristics at low 
and high lift coefficients. At moderate lift coefficients (0.6 to 1.0), 
it appears that LEEs with lower sweep angle are more effective in 
reducing the drag. The same conclusion can be derived from Fig. 3.6, 
where lift-to-drag ratio is examined at different lift coefficients. 

It appears instructive to compare the aerodynamic effectiveness of 
different LEEs relative to their planform area by considering the effect 
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of other geometrical parameters such as chord and sweep angle simultan- 
eously. In order to demonstrate these effects, a multi-variable plot 
shown in Fig. 3.7 was generated. This figure incorporates two aero- 
dynamic variables (i.e., L/D, a) and three geometrical parameters (i.e., 
sweep angle, constant chord, area) in a single plot. In general, the 
figure shows that the LEE planform area does not have a considerable 
effect on lift-to-drag ratio over the entire range of angle of attack. 
Further, with regard to the comparison of the aerodynamic effectiveness 
of LEEs with different constant chord and sweep angles, the following 
conclusions are drawn based on equal LEE planform area. 1) At moderate 
angles of attack (6° to 10°), it appears that constant chord LEEs pro- 
duce a better lift-to-drag ratio. 2) At 12° angle of attack, LEEs with 
sweep angles 57° to 55° generate better L/D, however outside this range 
constant chord LEEs achieve either the same or better improvements. 3) 
At 14° to 16° angle of attack, only low sweep angle LEEs appear to be 
more effective. However at higher angles of attack (18’ to 20°), the 
figure shows a very slight change in L/D ratio, regardless of the LEE 
planform shape or area. 

Evidently, these results coupled with the present study's design 
requirements (i.e., LEE with minimum area and chord) suggest that, 0.8" 
constant chord LEE has the aerodynamic potential of being selected as a 
candidate for the final LEE planform design. Further, 1.2" constant 
chord LEE also appears to be promising. Although this LEE benefits from 
1.5 times larger area compared to 0.8" constant chord LEE, it produces 
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14.3% improvements in L/D at 10° and a slight increase at 8 and 12 
angles of attack. 

3.4 Span Extent 

The selected LEE planforms which were used earlier in Secs. 3.2 and 
3.3 to investigate the aerodynamic effect of constant chord and sweep 
angle, are employed here to study the same effects on reduced span 
extent. The selected LEE span extents included 12.11" and 8.07" which 
correspond to 75% and 50% of the wing-model semi span, respectively. 
Associated with each span extent, twelve LEE planforms (i.e., 4.8", 

3.6", 2.4", 1.6", 1.2", and 0.8" constant chord, and 53°, 54°, 55°, 56°, 
57°, and 58° sweep angle) are examined in this section. The VLM-SA 
estimates were generated by using the same number of horseshoe vortices 
on the wing model, however, 7-by-18 and 7-by-12 array of horseshoe 
vortices were used on LEE planforms with 75% and 50% span extent, re- 
spectively. 

The VLM-SA analytical estimates of lift and drag coefficients for 
75% and 50% span-extent LEEs with constant chord are plotted in drag 
polar form in Figs. 3.8 and 3.9, respectively. These two figures coupl- 
ed with Fig. 3.2 (i.e., 89% span-extent LEE) show that in general beyond 

C , the reduction in the extent of the LEE's span causes an increase 
L *u 

in drag at constant lift coefficients. It is interesting to note that 
the drag increase associated with a longer chord LEE is relatively 
substantial. Similar drag polar plots are shown in Figs. 3.10 and 3.11, 
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Fig. 3.10 Effect of LEE sweep angle on drag poltir, n|FE=75% 
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Fig. 3.11 Effect of LEE sweep angle on drag polar, t^ee 


respectively, for 75% and 50% span-extent LEEs with different sweep 
angles. These figures coupled with Fig. 3.5 (i.e., 89% span-extent LEE) 
also show a similar drag increase behavior at constant lift coefficients 
as the span-extent of the LEEs are reduced. Further, these figures show 
that the drag increase associated with LEEs having lower sweep angle is 
relatively higher. The following paragraph examines the aerodynamic 
effectiveness of different LEE planforms as a function of geometrical 
parameters by including the effects of LEE span extent, sweep angle, and 
constant chord. 

Lift-to-drag ratio for 75% and 50% span LEE with different constant 
chord and sweep angle are plotted in Figs. 3.12 and 3.13, respectively. 

In general, a comparison of these figures with Fig. 3.7 shows that, a 
reduction in the LEE span extent decreases the L/D ratio in the a range 
of 6° to 12°, regardless of the LEEs planform shape and area. However, 
this effect appears to be insignificant as a increases beyond 14°. In 
regard to the comparison of the aerodynamic effectiveness of LEEs with 
different constant chord and sweep angle having the same planform area 
at a reduced span extent, the following conclusions are drawn. 1) 75% 
span LEEs (Fig. 3.12) conform with the conclusions drawn from 89% span 
LEE (see sec. 3.3, Fig. 3.7), and therefore they are not repeated here. 

2) At 6 8 to 10° angle of attack, Fig. 3.13 shows that 50% span extent 
LEEs with constant chord appear to be aerodynamical ly slightly more 
effective. However, in general, the L/D variation as a function of LEE 
planform shape or area appears to be negligible throughout the examined 
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Fig. 3.13 Effect of LEE-planform geometrical -parameters on 1 ift-to-drag 
ratio, n, ^,-=50%. 


a range. This is an interesting result. As shown in Fig. 3.13, 
although 3.6" constant chord LEE has 4.5 times larger planform area than 
the 0.8" constant chord, they both produce almost the same L/D ratio 
throughout the angle of attack range. Further, Fig. 3.14 shows the re- 
duction in L/D ratio as the LEE span extent is decreased, for two dif- 
ferent angles of attack, on the extreme LEE planform shapes. These re- 
sults suggest that a minimum LEE chord at the inboard of the wing lead- 
ing edges (i.e., apex region) is sufficient to yield an aerodynamical ly 
efficient LEE planform shape. 
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Chapter 4 


RECOMMENDATIONS AND CONCLUSIONS 

The present study demonstrated the applicability of a newly devel- 
oped analytical design procedure for the determination of an aerodynam- 
ical ly efficient Leading-Edge Extension (LEE) for thick delta wings. 
Although, the procedure is general enough to incorporate the effects of 
pitching moment, this study considered only the aerodynamic performance 
of a cambered and twisted wing at high-subsonic speed. Through an exam- 
ination of the available analytical tools, the PAN AIR and VLM-SA com- 
puter codes were employed to carry out the first application of the de- 
veloped design procedure for the given wing of the present study. The 
following sections simmarize the effectiveness of the analytical codes 
employed in this study as well as the resulting LEE effects. 

4.1 Analytical Tools Effectiveness 

1. The Pseudo-Stagnation Stream Surface (PSSS) solution associated 
with the wing model at the attached flow design condition was determined 
by employing the PAN AIR code. 

2. The PAN AIR code was further employed to reexamine the accuracy 
of the determined PSSS solution. This investigation showed that re- 
solving the problem with the PSSS in place resulted in a negligible 
influence on the aerodynamic pressures of the wing model at the attached 
flow design condition. Hence, it was concluded that the determined PSSS 
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solution is a good approximation of the dividing streamline surface. 

3. An evaluation of the analytical capability of VLM-SA demon- 
strated that, at low-to-moderate lift coefficient, the code estimates of 
drag were in good agreement with experimental data for a thick delta 
wing. However, when a LEE was added as shown in Appendix B, the VLM-SA 
code over-estimated the drag in the lift coefficient range of about 0.05 
to 0.80, because of the omission of suction pressures on the wing 
thickness. As a result, by analogy, it was concluded that the code 
estimates for drag would be too high in the present study of the wing- 
LEE configuration. Hence, throughout this study the aerodynanic effec- 
tiveness of different wing-LEE combinations were emphasized relative to 
one another rather than to the basic wing model. 

4.2 LEE Aerodynamic Effectiveness 

1. The analytical estimates obtained from VLM-SA code indicated 
that increasing the LEEs' constant chord reduces the total drag at mod- 
erate-to-high lift coefficient. However, no considerable improvement in 
L/D ratio was experienced with increasing constant chord size (i.e., in- 
creasing LEE pi anform-area) for the majority of the examined angles of 
attack. 

2. At low and high lift coefficient, the variation of the LEEs' 
sweep angle had only a slight effect on the drag and lift characteris- 
tics of the wing-LEE combination. However, at moderate lift coeffi- 
cients (i.e., 0.6 to 1.0), it was found that LEEs with lower sweep 
angle were more effective in reducing the drag. In general, based on 



equal planform area, it was evident that constant chord LEEs produce 
either the same or better L/D ratio than LEES with different sweep 
angles for most of the examined angles of incidence. 

3. Outboard reduction in the LEE span extent decreased the L/D 
ratio in the a range of 6.0* to 12°, regardless of the LEEs planform 
shape and area. However, this effect appeared to be insignificant 
beyond 14° angle of attack. 

4. The design requirements of the present study, coupled with the 
results obtained for thirty-six different LEE planforms, suggest that 
the 0.8" constant-chord with 89% span extent has the aerodynamic poten- 
tial of being selected as the best candidate for the final LEE planform 
design. Further, relative to the above LEE planform, the 1.2" constant- 
chord, with the same span extent, also appeared to be promising in 
achieving a better L/D ratio, especially at 10° angle of attack. 

4.3 Concluding Remarks 

It is firmly believed that the LEE concept has extensive drag-re- 
duction capability that justifies further investigation. The principal 
research effort should be toward the selection of appropriate analytical 
tool(s) for the accurate determination of the PSSS at a given design 
flow-condition, as well as, the LEE planform optimization. Although 
efforts have been made to validate the results obtained in the present 
study wherever possible, it is recommended that the aerodynamic 
effectiveness of the resulting LEE be verified experimentally in a wind 
tunnel investigation. 
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APPENDIX A. PAN AIR CODE 


The intended purpose of this Appendix is to discuss briefly some of 
the capabilities and the limitations of the PAN AIR code. Also, the 
degree of accuracy of the determined PSSS solution is examined. 

A.l Discussion 

PAN AIR code is a boundary-value problem solver for the Prandtl- 
Glauert equation 


( 1 -M 2 )<J> + <f> + <f> = 0 

Thus, PAN AIR's nunerical solutions embody all the limitations inherent 
in this approximate representation of the flow field. Such limitations 
include steady, inviscid, irrotational , subsonic and supersonic flows. 
This higher order panel code uses linear source and quadratic doublet 
strength distribution to determine the flow field. 

PAN AIR aerodynamic analysis capability consists of ability to: a) 

calculate pressure distribution and velocity components at any point on 
the surface of a configuration and b) calculate forces and moments both 
on the configuration as a whole and on specified portions of the config- 
uration. This potential flow code is capable of employing survey net- 
works which can be used to determine the velocity components at any 
point away from the configuration surface. In addition to the above, 
the PAN AIR code has the Influence Coefficient (IC)-update capability. 
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This capability, enables one to examine configurations which differ from 
the ones already processed in a limited fashion with respect to geom- 
etry. In this process a configuration surface is partitioned into 
several networks, with one or more tagged "updatable" in the original 
submission. The subsequent run, with any change (i.e., size, location) 
in the updatable networks, can utilize some of the unaffected 
calculations which have been performed and saved from the original 
execution. 


A. 2 PSSS Validation Test 

As discussed in Chap. 2 of the present study, the PSSS solution was 
determined for the wing model at the attached flow design condition by 
the PAN AIR code. As part of the present study, it appeared necessary 
to examine the accuracy of the determined PSSS solution. This test was 
conducted by employing the PAN AIR code once again to model the wing- 
PSSS combination at the same attached flow design condition used in the 

S 

present study to determine the PSSS solution. Although the modeled PSSS 
is specified as being a lifting surface. Fig. A.l shows that its pres- 
ence has little effect on the pressure distribution over several typical 
sections of the wing upper and lower surfaces. Also, as shown in Fig. 

A. 2, the lifting pressure across the PSSS appears small at the same 
typical sections, especially away from the PSSS leading-edge. From 
these results, it is evident that the attachment of the PSSS solution, 
does not have a considerable influence on the performance of the wing 
model at the attached flow design condition. Therefore, it is concluded 



that the determined PSSS solution is close to the actual dividing stream 
surface (i.e., pseudo-stagnation stream surface). 
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APPENDIX B. VLM-SA CODE 


This Appendix includes a discussion on the VLM-SA methodology by 
highlighting some of the code capabilities and limitations in Sec. 

B.l. Also, the computer programs developed to generate the required 
local angles of attack for twisted and canbered wing and the LEE surface 
are listed in Sec. B.2. Lastly, results and a discussion on VLM-SA 
evaluation is presented in Sec. B.3. 

B.l Discussion 

Theoretically, the conventional vortex lattice method is incapable 
of determining the aerodynamic characteristics of a wing configuration 
having leading-edge vortex flow. However, to overcome this limitation 
Lamar and Gloss [15 j of NASA Langley have developed a computer code 
which couples the conventional method with suction analogy concept (VLM- 
SA). The suction analogy [17] states that the attached flow leading 
edge suction force which no longer acts in the chord plane when the 
leading-edge flow separates, is reoriented normal (rotated 90°) to the 
upper surface by the vortex flow action, thereby producing an additional 
force. This force is the aerodynamic force associated with the leading 
edge vortex flow. The VLM-SA code calculates this force and its contri- 
bution to lift, drag, and pitching moment. Then, these contributions 
are added to the potential flow results to find the total vortex induced 
forces and moments. 



The VLM-SA code was extended to account for twist and camber (call- 
ed VLM mark 4.0 version) as reported by Lamar/Herbert in refs. 18 and 
19. VLM mark 4.9, the latest version of VLM-SA code, is capable of 
accounting for the effect of leading-edge radii as well. However, it 
does not provide pressure loading information. Although, the code 
neglects the effect of thickness, there have been techniques devised to 
include them [20]. 

The latest version of VLM-SA code is employed in the present study. 
The code aerodynamic analysis capability consists of ability to estimate 
overall potential flow, as well as, vortex flow induced forces and mo- 
ments on complex planforms at subsonic speeds. Figure B.l, shows three 
types of flow situations considered by the code. These are: 1) full 

leading-edge suction (attached flow), 2) zero leading-edge suction 
(attached flow), and 3) zero leading-edge suction (vortex flow). All 
predicted results presented in this thesis correspond to the zero lead- 
ing-edge suction with vortex flow. 



Fig. B.l Theoretical flow mechanisms employed in the extended VLM-SA 
computations. 
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•Hit'S LOCAL CHORD. 

KLB>2SR 057 CHORD STATIONS. 

-FRACTION OF WINS LOCAL Q-®\9. 

-WING SEMISPAN. 

-PRACTICE OF tflhQ SEHIS^AN. 

-WING ROOT CHORD. 

5 t24r ( lj£fi?i£; v S-OPSS (CHORDWISE HORSE- 

T0 65 AT A SPAN STATION. 

1 ' DESIRED NUMBER OF SPAN STATIONS AT WHICH 
HORSESHOE VORTICES WILL BE LOCATED. 

J N ^-COORDINATES FOR SLOPE CALCULATIONS. 
-^-COCKDINATE ALONG THE MEAN CAMBER LINE JUST BEFORE EACH CONTROL 

-^-COORDINATE ALONG THE MEAN CAMBER LINE JUST AFTER EACH CONTROL 


u^rwwVre 11 ^ ?S?§ h lII?l A S SS $£i A S WITH THE 'LIBRARY SUBROUTINE < ST I IN I ) 

THIS PROGRAM. The FIRST, INTERPOLATES THE SLOPES AT THE 
CONTTOL POINT OF EACH GRIDE LINE. THE SECOND, INTERPOLATES THE SLOPES 
a^CWTROL POINT OF EACH ELEMENTAL PANEL FROM THE SLOPES CBTAIhCD 

h!' tS ‘iWX INTEGER SPECIFYING THE NUMBER OF NODES. 

Ml, M2 - INPUT INTEGER SPECIFYING THE NUMSER OF STATIONS TO BE INTERPOLATED. 
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OOUUOOUUUUUOOOUUOUUUUOOOOOOOOUOOUOOOUOUU 


Ml. M2 
XLWI 

ZUNI 

TUNI 

SIGMA 


lEtcat 

END 

IV 

FUNI 


VKU 

IERR 

XCPT 

VCPT 

L SLOPE 

FSLQP 


SLSLOP 

SLOPCP 

SLOPES 


ALPHAL 


■INPUT INTEGER SPECIFYING THE NUM3SR OF STATIONS TO BE INTERPOLATED. 
■INPUT REAL ARRAY OF LEN3TH N1 CONTAINING Tf€ X-COORDII^TE OF 
THE NODES. < NONDIKENSICfGIZED WITH LOCAL CHORD.) 

array of length ni containing the z -coordinate of 

THE NODES. (NONDIKENSIGNLIZED WITH LOCAL CHORD 7)^ 

-If^REAL ARRAY OF LENGTH Ml CONTAINING n«XH*OORDINATE 
t£ INTERPOLATED VALUES OF 

TNG Z -COORDINATE ARE DESIRED. 

-IT IS A REAL NUPSER, CONTAINING THE TENSION FACTOR. SIGMA 
INDIMTES THE CUR VI NESS DESIRED FOR INTERTOLATION? ZERQSIGHA 
PRODUCES AN EXACT CU3IC S*LUE FIT. WHILE IF SIGMA IS FIFTY 
TNG RESISTING CURVE IS PIECEWISE LINEAR. 

■JtfVT INTEGER SPECIFYING THE C ONDITION ON THE CURVE AT THE 
TWO EM) POINTS. IEM>SW-2 CORRESPOM)S TO Y“«070 . 

"OUW VIABLE. NOT USED IN THIS PROGRAM. 

-INPUT/OUTPlfT INTEGER. 

O UTPUT A RRAY OF LENGTH Ml CONTAINING TVE VMJJES OF TNC 
INI^^ATg Z-COg^IN ^TES CORRESPONDING TO TUNIX. ( I.E. TUNIX. 
S£?Ig® ££*A AW Z COO RDINATE S OF TKE CONTRO L POIN TS ALONG 
TIG KEAN CAMBER LINE OF DIFFERENT WING SCTION, RESPECTIVELY. ) 
-INPUT/OUTPUT WORK REAL ARRAY OF LENGTH 2N1-1. 

^^TPUr 1WEGSR .ERROR CODE. ZERO IERR CORRESPOWS TO NORMAL RETURN. 
-LO CAL S LOPES IN CHORDWIS DIRECTION. 

-Y-COORDIKATE OF THE CONTROL POINT LOCATIONS MJONB TVC MEAN 
CAHEER S URFAC E AT WHICH THE SLOPES ARE DESIRQL^ 

■Ugjffi-, ^ PIRSCTION alow tvg gride lines, starting 

FROM 10 TE PROCEEDINS from iwqard to outbomd. 

*c^^c.%SSJ N «SSf® tt i^S IfSC I12! AT control point of 
E^HD-BGNIT^P^L, STATING from UCOAflD TO OUTBOARD AM) 
PROCEEDING from LE To te. 

"15^.Sy>Ss IN CHORDWISE DIRECTION ALONG T>€ GRIDE LUCS. 

S^TIW FRCHIN£OAg>JU OUTBOARD AM) PROCEEDING FROM LE to TE. 

®T HE SWZ AS ABOVE FSLGP. 

^9££kJy® s I?.JN 5Pft?vWI SE _ D I^CT I ONAT THE CONTROL POINT OF EACH 
Ej^^^Pg^jrART 1 !® FROM LE TO TE AM) PROCEEDINS FROM 

-LOCAL ANGLES OF ATTACK (SLOPES) AT THE CONTROL POINTS IN RADIAN. 


P^S^.LALPI^ INPUT, OUTPUT, TAPES, TAPES, TAPE! 0 ) 

« Z cuh§? > £3£tf! I6,G0),LCHORD(18,M),CHDfiD< l6, 1 >,PSPWJC 16) 
),TLNIX<14), W<U< 1 10), XUNKG0), YUNKS3), 
i 6 i i 4 } ' PUNIP’ <16.14). WXU2( 31 ) , SL0PES<64, 14), 

ipiloSfji?: hi. 

INTEaR ) IE»C)S»<2>.I«XS.0rS 

R2K !^^'!^^<a 5& T/2S * 33709/ ' J1ST/1/, SPAN/18. 147G0/ 
g*TA 2 1 /’ytIg£SW / 2, 2/, SIGMft/33;/,DNXS/14/,DYS/24/ 

g£TA S^/a"S^i'S*§9S6^ ,6/ '^ 24/ t2^N 0/,4 ‘ / 

*1™' ‘HISS' -§§122' -40300, • 448S3 ' .480®. .52500, .57000, 
-flggg- -68080. .733C0, .77853, .82727, .86727, .00363, 

I • Ww30 , 4 9 && 0 Q/ 

DO 10 I -1,NS 


69 



le 


11 


12 


13 


17 

IS 


111 

lie 

112 

120 

160 


266 

14 


161 


DO !G I »1,NS 
DO 10 J °1 , NODE 

6 ?* TC 1 Its* > «■ * T < I * J i 

word I. > -|Ro6T-xru,n 
DO M I-1.NS 
DO II J»I,NOO£ 

HQKH1 “NODE"! 

DO 13 I»1.16 

2CI. 1 >«ZBt 1,1 > 

ZNCAM 1.1 )•£(!, n 
LCKORDti,!)-©.' 

DO 12 J-2.N0DEMI 
JH1»J“1 

« 6 do« .. , , 

LCHCE2X I,N0DE>«1. 

Z( I t N0DS)»2B< I JfflOE) 
agfejn^WBE )«$B< I , NODE ) 

IttO 

DO 14 I«*1.N3 

DO 1GD K=1 ,2 
DO IS -M.KOBS 
IF(K.Efl.2) 60 TO 17 
»MIU)«lD{Qr3>U,J> 

YUNH J)«B££Kt I. J)/OJQf®( 1 . 1 ) 

IFU.GT.DNXSIOO TO 15 

ra5Ti5ail” <! FLOft1 ' 1 J: ‘-s* 25 >/chotono*constck) 
i iej«> st Iuru<N, ' ni ' x,J?,I ‘ YIWI • tunix, siem, id gov, em>, iv, funi , t*u, 

IF( IERR.KE.0) GO TO 21 
IFOC.EQ.DGO TO 111 

IF<K.£Q.2) GO TO 112 

w lie J=I,DNXS 

CONTINUE 
GO TO 1C3 
M 120 J»1,DNXS 
RWIP1(I.J)-RJNI<j) 

CONTINUE 

CONTINUE 

DO 2^ J»1,DNXS 

Sn^l' J >SB( RWIP1 ( 1 ' J , " RJNIM1 * 1 ' J ’ ^ 082 ’ 

CONTINUE 

DO 161 Jal^DNXS 
DO 161 I-I.NS 

IW*6 

00 102 J-1.DNXS 



103 


106 

102 


C 

C 


2 

21 


DO 102 J-1.DNXS 
DO 103 1-1. NS 
xgr;i ii |SLa.fip<j.:, 

1F^ S ^M^^' PSf ^' XCPT ' VCPT ' SI0HA ' lETOSV'B®. IW, 

tF< ieJr.ic!©) GO TO 21 
DO 105 I-I.OYS 
a^PCPU. I^-FSLOP(I) 

CONTINUE 
CONTINUE 
DO 300 I-1.DYS 
DO 300 >K0NXS 
SLOPESI I , J )«SL0PCPU, I) 

J>«- atancsl6pesi i. JM 

"* CeTOT,I * D SL0PE3 

DO 600 1*1 ,OYS 
K*OYS*l-I 

Sntxnue' 2,< ^ lphw ’ ck * Jl * J-! « 0NXS > 

FORMAT* 7F1O.0) 

REWIND 10 
STOP 
EM> 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


i r ssssss? 


; RATE &\LY THE LOCAL AN&.ES 

* GENTRIES AS REaJIREDBV T>€ VLM-SA Si m 

* IS SLMMERIZED IN THEFCLLOyi^ TWOSTEPS ' ™ IS PR0QRAM 

* j ] DETERMINING T(£ SLOPES ALONG EACH SPAN STATION 

* i^jLjyTE^OLATING TO DETERMINE THE DESIRED NUMBER OF 

* nif^ A suSACE? NTR0L PI0NT ** EACH el£ment ^^nelon * 

vie v 7i £ -w 


XLE. Y, 
NONZLE 
LECORD 
LECORO 
PL SPAN 


— LIST Or VARIA&£S 

Z^lfcpTE^fcm. CHCKD. 
5rcS r §«S. L0CAL CH0R0 * <M*®IK3«I0?fi.IZED WITH LOCAL CHOf©) 

^CTIg^g > LOCA L SPAN. <N0M)IK SM3I0MLIZED WITH WING SEHISFAN) 


10 


11 

12 


15 
17 

16 


««< CtKSfTS FOR THE ISmiNge^l^SRE GIVEN IN THE PREVIOUS »>» 
PROGRAM LJEESLP ( INPUT, OUTPUT, TAf^S TAPEIfl) 

i^Ronl! i^lleco^j !?j' wrCLEt M « ,3) ' 

INTEGER i£rD3W< 2 ) , 6nXS. DYS 

^]y.NSDHyia/.S>AN/16. 147ES/ 

as 

! -§I^' 

DO 10 >1 'f33>£ 

ro^i 1 ' *i 2 i Its' J ’ • Y< 1 • J 1 1 2LE ‘ 1 • J » 

PLSPAN(I)-Y(I, 1 1/SPAN 
DO 11 J=! .NODE 

LEC^J I )»xl£( I . NODE )-XLE( 1 , 1) 

J,-XL£< I. 1 > >/LECGRDt I ) 

N»JZLEU, J)«vZLE(I, J)/LECQHD(I) 

CONTINUE 
IW=0 

DO 14 I»|,ns 
DO 108 K*i; 2 
DO 15 J*1,NOOE 
IRK.EQ.2) GO TO 17 
XUNI< J)»PLCQRD( I, J) 

^UJWBNZU^I... 

DO 16 M»J,7 

cSItinue"' FL 0AT (M, " 0 - 25 ,/( CHORDNO > -CONST C k ) 

CALL STIUNKN1 , Ml ,XUNI,YUNI. TUNIX, SIGMA. IENOSW.EM), IW,RJNI,WKU, 


1 10 
112 

120 

ta 


14 


600 

2 

21 


I JEW> STlLWI(N, ' M1 ' XUNX#riJNI ' TUNIX#SXGHA# IENDStf,EN D , IW.FUNl.ttU. 

IF(IERR.NE.O) 00 TO 21 
IFIK.EQ.I >G0 TO 111 
... i£«K.EQ-2> 60 TO 112 
Ml 00 118 J-I.DNXS 
fwihhi.jj-Aunhj) 

CONTINUE 
GO TO 100 
00 (20 J-I.DNXS 
FJMPIU.JJ-^NHJ, 

CONTINUE 
CONTINUE 

00 200 J-I.DNXS 

LgjOPEn. J)«<FIWIP1 ( I. JI-FUNIHI ( I. J) )/<2*.0O2) 

CONTINUE 

00 101 J-I.DNXS 
00 101 I-I.NS 
CLSLOPU. I)-LSLOPE(I.J) 

101 CONTINUE 

iw-e 

00 102 J-I.ONXS 
DO 103 1-1. NS 

xcpt(I)<lsl6pu.i> 

CONTINUE 

1 rao> S iaS' i&T *" PLSfm - xa ‘ T ■ ' nr * ■ si—. IOCS*. ew>. Ik. 

ffaftSift* ro 21 

CONTINUE 
DO 300 1-1. OYS 
DO 300 >1,0NXS 
SLgPESU,J)-SLOPCP(J t I) 

J>»- ATAN< Sl6pES< I.J)) 

00 600 I-1.DYS 
K«0YS*1-I 

«ITEn0.2)(ALPHflL<K, J). J-1 .DNXSl 

FORMAT( 7 F 10 . 4 ) 

REWIND 10 
STOP 
EM) 


103 


106 

102 
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B.3 Evaluation 


As a part of the present study, it was important to examine the 
analytical capability of the VLM-SA code for thick wing configuration 
with leading-edge extensions. For this purpose, the experimental data 
obtained by Wilson and Lovell, on the thick DM-1 with and without 
leading edge extension, was selected for validating the results obtained 
from the VLM-SA code. Although the effect of leading-edge radii is in- 
cluded in the resulting VLM-SA solutions, the thick DM-1, which is a 
synmetrical wing configuration with 0015-64 NACA airfoil section and no 
twist, is approximated by its projected planform (flat DM-1) in this 
analytical study. Experimental lift data obtained by Wilson and Lovell 
on the DM-1 with and without the leading-edge extension, as well as, the 
resulting VLM-SA for the same configurations are shown in Fig. B.2. 
Obviously, the code solutions over-estimate the lift for both the DM-1 
and DM-1 + LEE combination throughout the angle of attack range. How- 
ever, as shown in Fig. B.3, the drag polar comparison shows that, in 
case of the DM-1, the VLM-SA solutions agree well with experimental data 
up to lift coefficient of about 0.6. Beyond this lift coefficient, the 
experimental data tends to deviate from the code's solution, because, as 
was reported by Wilson and Lovell, the flow over the basic DM-1 appears 
to become disorganized and result in an increase in drag and a decrease 
in lift. As a result, the drag polar curve associated with experimental 
data is higher than the estimated VLM-SA solution. In the case of the 
DM-1 + LEE combination, the VLM-SA over estimates the drag in the lift 
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Fig. B.2 Theoretical and experimental! ift characteristics. 
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Fig. B.3 Theoretical and experimental drag polars. 


coefficient range of about 0.05 to 0.80. This difference was rather 
expected, because the resulting VLM-SA solutions do not include the 
effect of the low pressures acting between the LEE and upper surface 
maximum thickness line of the wing section to produce a thrust. Hence, 
the computed Cq values are higher than the experimental data. This 
effect can be seen in Fig. B.4, where lift-to-drag ratio is plotted 
against lift coefficient. This figure shows again a fair agreement 
between the code solution and the experimental data up to lift coeffi- 
cient of about 0.6, for the basic DM-1. However, the L/D theoretical 
curve associated with the DM-1 + LEE combination is considerably lower 
than the data points, because of higher drag estimation (note that the 
lift is also over estimated by the code. Fig. B.2) by the code. There- 
fore, by analogy it is expected that the VLM-SA solution of drag would 
be higher in the wing-LEE analysis for the present study. 



Fig. B.4 Theoretical and experimental lift-to-drag ratios 
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